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ABSTRACT 

Gamma Ray Bursts (GRBs) show evidence of different light curves, duration, afterglows, host galax- 
ies and they explode within a wide redshift range. However, their spectral energy distributions (SEDs) 
appear to be very similar showing a curved shape. Band et al. (1993) proposed a phenomenological 
description of the integrated spectral shape for the GRB prompt emission, the so called Band func- 
tion. In this letter we suggest an alternative scenario to explain the curved shape of GRB SEDs: the 
log-parabolic model. In comparison with the Band spectral shape our model is statistically favored be- 
cause^ it fits the GRB spectra with one parameter less than the Band function and it is motivated by a 
tlicoretic;aI acceleration scenario. The new Fermi observations of GRBs will be crucial to disentangle 
between these two models. 

Subject headings: stars: gamma-ray burst: general, radiation mechanisms: nonthermal, acceleration 
of particles. 



1. INTRODUCTION 

The physical mechanisms behind the GRBs prompt 
emission are still under debate. Band et al. (1993), 
investigating the BATSE GRBs sample proposed a 
phenomenological description of the integrated spectral 
shape for the GRB prompt emission, the so called Band 
Junction. The introduction of this function was strongly 
suggested by the observational evidence that the shape 
of the Spectral energy Distribution (SED) of the GRB 
prompt emission is convex and broadly peaked. It is re- 
markable that there has not been physical explanation 
in terms of accelerations processes and non-thermal ra- 
diative losses that can led to the Band spectral shape. 

In this letter we propose to describe the curved shape 
of GRB prompt emission using the log-parabolic function 
(Massaro et al. 2004), successfully used to describe the 
SEDs of BL Lac objects over several decades. First, we 
consider the differences between this model in compar- 
ison with the Band function investigating the different 
7-ray flux predictions in the Fermi LAT energy range. 
Second, we point toward the physical interpretation of 
the log-parabolic shape in terms of Fermi acceleration 
mechanisms. Finally, we present a simple synchrotron 
emission model to explain the GRB prompt emission, 
that appears to be the most reasonable scenario. 

For our numerical results, we use cgs units unless 
stated otherwise and we assume a flat cosmology with 
Ho = 72 km s-i Mpc-\ = 0.27 and I^a = 0.73 
(Spergel et al. 2007). 

2. THE SHAPE OF THE SPECTRAL ENERGY 
DISTRIBUTION 

GRBs have a non-thermal spectrum that varies 
strongly from one burst to another. It is generally found 
that a simple power law does not fit well their spec- 
tra because of a steepening toward the high energies. 
The Band function phenomenological model (Band et al. 
1993) describes the prompt time- integrated GRBs spec- 
tra, composed by two power laws joined smoothly at a 
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where F{E) is the number of photons per unit of area and 
energy and time, while Eq is a reference energy usually 
fixed to the value of 100 keV. Under the continuity re- 
quirement for the function F{E) and its first derivative, 
the break energy and normalization are given by: 
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There are no particular theoretical scenarios that pre- 
dict this spectral shape making it only a phenomenolog- 
ical model. However, it provides good fits to most of the 
observed spectra in terms of four parameters, namely: 
the two photon indices a and j5, the exponential cut-off 
Ec and the normalization constant Fq, with all four pa- 
rameters directly estimated during the fitting procedure. 
The peak energy Ep of the SED (i.e. S(E)=E2F(E)) is 
related to the spectral parameters by: 

(4) 



Ep = {a + 2)E^ < E^ 



for which, typically, — 2 < a < — 1. We also note that 
for typical values a ~ —1.4 and /3 ~ —2.4 (Band et al. 
1993), Eb « E^. 

We propose to describe and interpret the shape of the 
SED in the GRB prompt emission using a model defined 
by the following equation: 



F{E) = Fo 
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where Fq is the normalization, a is the spectral index 
at energy Eq and b is the parameter which measures the 
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spectral curvature. This model is known as log-parabolic, 
a curve having a parabolic shape in a log-log plot (Mas- 
saro et al. 2004, Massaro et al. 2006). We remark that 
this spectral distribution is the classical log-normal sta- 
tistical distribution. In particular, for this function, it 
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Fig. 1.— Two SEDs of GRB 910601 (upper panel) and GRB 
920622 (lower panel) with the log-parabolic best fit model (black 
line). It is clear how the spectral shape is well described by the 
log-parabolic function (data taken from Schaefer et al. 1994 and 
from Tavani et al. 1996). 

is possible to define an energy dependent photon index 
T{E) given by the log-derivative of Eq. (5), 

T{E)^a + 2blog{E/Eo) (6) 

which describes the continuous change in the spectral 
slope. 

The peak energy Ep and the height of the SED S(E) 
calculated at its peak frequency 5*^ can be evaluated by 
the following relations 

Ep = EolQ^ (7) 

5p = 5*010^ (8) 

where 5*0 is S{Eq). 

Consequently the spectral shape can be expressed in 
terms of 6, Ep and Sp using the relation: 

S{E) = E^F{E) = Sp 10-'' i°g'(^/-^=p) , (9) 

where Sp — EpF{Ep). In this form the values of the 
parameters 6, Ep and Sp are estimated independently in 
the fitting procedure, whereas those derived from Eq. (7) 
and Eq. (8) are affected by intrinsic correlations (Tani- 
hata et al. 2004, Tramacere et al. 2007). 

As example of the goodness of the fitting procedure, we 
show in Fig. 1 the best fits of two GRB SEDs evaluated 
with the log-parabolic model namely GRB 910601 and 
GRB 920622, two of the most bright and well studied 
GRBs present in literature, (data taken from Schaefer 
et al. 1994 and Tavani et al. 1996). Their best fit pa- 
rameters are: Ep = 0.47 ± 0.01 MeV, b = 0.74 ± 0.03, 



Sp = 9.99 ± 0.12 X 10-*^ erg s^^ cm^^ for GRB 910601 
while Ep = 0.460 ± 0.003 MeV, b = 0.95 ± 0.02, Sp = 
1.098 ± 0.005 X 10-^ erg s'^ cm'^ for GRB 920622. 

We note that GRB spectra appear to be narrower with 
respect to the BL Lac objects (Massaro et al. 2008), 
having the curvature parameters close to 1. 

3. LOG-PARABOLIC VS BAND MODEL 

From a statistical point of view, the log-parabolic 
shape requires one parameter less than the Band func- 
tion and so is favored. There are two main differences be- 
tween these two models: first, the slope at low energies of 
the Band function is a power-law while the log-parabolic 
one has a milder curvature and second, the high energy 
tail is naturally curved and the expected flux in the 7-ray 
band is lower than the one predicted by the Band spectral 
shape (see Fig. 2). In particular, the log-parabolic model 
can describe a continuous curvature over the whole spec- 
trum while the Band function can only mimic it around 
the SED peak. 





• GRB 920622 




— Fermi LAT sensitivity 




— log-parabolic spectrum 




— Band spectrum 




■ ■ ■ Band spectrum + high energy cut-off 































10"' 10" 10' 10" 10"' 10* 

E (MeV) 

Fig. 2.— The SED of GRB 920622 and the comparison between 
the Band and the log-parabolic function used to describe its emis- 
sion. The extrapolation of the Band function predicts that a simi- 
lar GRBs should be detected within 100s in the Fermi LAT energy 
range, while the log-parabolic model is more in agreement with the 
GRB detection rate of the first year of Fermi LAT observations. 
The Fermi LAT sensitivity has been evaluated from that reported 
in Atwood et al. (2009). 

In the recent Fermi observations only 9 GRBs have 
been detected at high energies in the LAT energy range 
(Granot et al. 2009) in comparison with the predic- 
tions provided by the extrapolation of the Band function 
(Omodei et al. '2007, Band et al. 2009, Omodei et al. 
2009). Several explanations have been proposed to cor- 
rect the expectations and a high energy cut-off has been 
introduced in the Band function to arrange the lack of 
the observed GRBs in the Fermi LAT band (Band et al. 
2009). The introduction of this exponential cut-off in- 
creases the number of parameters in the Band function 
while the log-parabolic model appears to have a natu- 
ral explanation for the Fermi observations without the 
introduction of any new spectral parameter. 

In Fig. 2 we plot the SED for the GRB 920622 with 
different spectral model extrapolations in the Fermi LAT 
energy range. The Fermi LAT sensitivity evaluated for 
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an exposure of 100 s and with three different backgrounds 
as reported in Atwood et al. (2009) is also shown. Using 
the Band model a typical GRB with a similar spectral 
shape and total fluence of GRB 920622, as those reported 
in the BATSE bright GRB catalog (Schaefer et al. 1994), 
is expected to be detected in the Fermi LAT energy range 
while the predictions of the log-parabolic shape are very 
different and no LAT detection is expected for this GRB. 

We fitted with the log-parabolic model the SEDs of all 
brightest GRBs present in the BATSE catalog (Schae- 
fer et al. 1994) detected during 1991. The shape of the 
SED is well described in terms of this model as already 
shown in Fig. 1. We found that extrapolating the log- 
parabolic spectrum only 2 GRBs, out of 20 total, are 
expected to be marginally detected in the Fermi LAT 
band within 100s, more in agreement with the detection 
rate of the first year Fermi observations with respect to 
the high Fermi LAT detection rate estimated using the 
Band model extrapolation (Omodei et al. 2009, Band et 
al. 2009). In Fig. 3, It is evident how the curved shape 
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Fig. 3. — The log-parabolic models of the 20 brightest GRBs in 
the BATSE catalog detected during 1991 (Schaefer et al. 1994) It 
is clear how the extrapolation of the log-parabolic function led to 
conclude that only few GRBs are expected to have a detection in 
the Fermi LAT energy range in contrast with the expectations of 
the Band model. 

described by a log-parabolic function is more in agree- 
ment with the Fermi LAT first year GRB detections. As 
in Fig. 2, the Fermi LAT sensitivity has been evaluated 
from that reported in Atwood et al. (2009) rescaled for 
100s. 

4. ACCELERATION MECHANISMS AND SYNCHROTRON 
RADIATION IN THE PROMPT EMISSION 

The energy spectrum of accelerated particles by some 
statistical mechanism, such as those occurring in shock 
waves, is usually written as a power law. The origin of 
this interpretation resides in the first order Fermi accel- 
eration mechanism (Bell 1978, Blandford & Eichler 1987 
and Protheroe 2004), originally presented to explain the 
cosmic ray spectrum. However, the observational evi- 
dence that the SED of BL Lacs objects has a curved 
shape has demanded toward a different interpretation. 
In particular. Landau et al. (1986) provided a useful de- 
scription of the synchrotron component of BL Lac objects 



in terms of a log-parabolic model, that has been recently 
applied to describe the synchrotron X-ray spectra of TeV 
BL Lacs (Massaro et al. 2008). 

The theoretical interpretation of the log-parabolic 
model resides in the energy distribution of the emit- 
ting particles. The general solution of the kinetic equa- 
tion, for the time-dependent distributions with respect 
to the energy, yields curved particle energy distributions 
(PEDs) in the form of log-normal function, where terms 
taking into account of the stochastic and systematic ac- 
celeration by the Fermi mechanisms are considered (Kar- 
dashev 1962). 

Assuming a simple S function as initial condition for 
the PED, the analytical solution of the particle kinetic 
equation yields the log-parabolic shape in the form: 

/ \ -s-r log (7/70) 

N{^)^Noi^] (10) 

where the parameters s, r and the normalization Nq are 
directly linked to the physical parameters Ai and A2, 
while 7o is a reference energy. We note that the PED 
curvature r is only directly linked to the diffusion coef- 
ficient, this means that the curved shape of the parti- 
cle distribution depends on considering the second order 
Fermi acceleration mechanisms (Kardashev 1962, Paggi 
et al. 2009). The case r = corresponds to the simple 
power-law electron distribution as expected by a first or- 
der acceleration mechanism. The mean quadratic energy 
of the PED, (7^), corresponds to the second normalized 
momentum and can be expressed in the form: 

(7^)=7olO(^-^)/^^ = 7. (11) 

which corresponds to the peak of 7^ By freezing 

the PED slope, s, to the value of 2, as expected in the 
first order Fermi acceleration mechanism, it is possible 
to describe its shape in terms of the PED curvature r, 
7p and the normalization A^o and Eq. (8) can be written 
in the form 

/ \ -2-rlog(7/7p) 

N{^)^Noi^^j . (12) 

Finally, we remark that it has been recently shown 
that including synchrotron and inverse Compton radia- 
tive losses as well as of the "disappearance" of fast parti- 
cles that escape from the acceleration region, either as a 
result of nuclear collisions or escape from the acceleration 
region, the numerical solution for the kinetic equation 
can be successfully described in terms of a log-parabolic 
shape (Tramacere et al. 2009, Paggi et al. 2009). 

Under this assumption that the synchrotron radiation 
is emitted by a log-parabolic PED (Eq. 10), the resulting 
flux density and consequently the SED can be well ap- 
proximated by a log-parabolic shape, expressed as: Eq. 
(5). An alternative scenario is based on the assump- 
tion that the probabilty to accelerate particles depends 
on energy in a simple relation given by P cx 7"'', and in 
this case the resulting PED yields toward a log-parabolic 
shape (Massaro et al. 2004, 2006). 

Applying the numerical code developed by Massaro et 
al. (2007) and presented in Paggi et al. (2009), we calcu- 
lated the synchrotron emission by a log-parabolic PED 
to describe the observed SED of GRB 910601. This is a 



4 



F. Massaro, J. E. Grindlay, A. Paggi, 2009 




E (MeV) 



Fig. 4. — The synchrotron model applied to the case of GRB 
910601. The intrinsic source parameters are reported in Table 1. 



TABLE 1 

The intrinsic source parameters for GRB 910601. 



Parameter 


Symbol 


units 


value 


rodshift 


z 




1.0 


FED slope 


s 




2 


FED curvature 


r 




7.58 


FED energy peak 


7p 




1.97X 10* 


FED minimum energy 


Tmin 




10^ 


FED maximum energy 


^max 




5x 10* 


electron density 


ne; 


cm 


1.13X 10* 


beaming factor 


5 




30 


magnetic field 


B 


Gauss 


10* 


volume 


V 


cm^ 


10*2 



clear example of how the log-parabolic scenario success- 
fully describes the shape of the GRB prompt emission. 
The parameters assumed for our calculations are given 
in Tab. 1, and in Fig. 4 we show the GRB 910601 SED 
with the model adopted. We fixed the redshift of this 
GRB to 1 because it is unknown. The parameters in 
Tab. 1 are all consistent with plausible values of the 
GRBs emitting region (e.g. Meszaros 2002 and zhand & 
Meszaros 2002). The synchrotron model, evaluated with 
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a log-parabolic FED, is in agreement with the data, ap- 
pearing to be a good description of the GRB SEDs of the 
prompt emission. 

5. GONCLUSIONS 

In the present letter we propose to interpret the SEDs 
of GRBs prompt emission using the log-parabolic shape. 
In comparison with the Band function (Band et al. 1993) 
the log-parabolic shape is favored for two main reasons. 

First, it is statistically better, because it requires only 
3 parameters, namely the curvature h, the peak energy 
Ep and the height of the SED evaluated at the peak 
energy Sp (see Eq. 6 and Eq. 15), while the usual Band 
model needs 4 spectral parameters, namely: a, /3, Ec 
and Fo or 5 if another high energy exponential cut off is 
introduced. Second, the proposed function has a strong 
physical motivation. This shape is directly related to the 
solution of kinetic equation for the particles accelerated 
by Fermi mechanisms when the random acceleration is 
also taken into account with all the other terms (e.g. first 
order Fermi mechanisms) (Kardashev 1962). 

In the recent Fermi LAT observations only few GRBs 
have been detected in contrast with the predictions of 
the Band function. The high energy curvature of the 
log-parabolic shape has a natural explanation for the 
Fermi observations without the introduction of any new 
parameter, as the exponential cut-off in the Band func- 
tion. 

As shown in Massaro et al. (2006) or more recenly in 
Tramacere et al. (2009) and Paggi et al. (2009), the 
synchrotron emission of a log-parabolic electron distri- 
bution yields a curved SED near its peak, well described 
in terms of the same spectral shape. In this letter, we 
also presented how the synchrotron scenario from a log- 
parabolic FED can describe the spectrum of the GRBs 
prompt emission, successfully. 

Finally, we remark that a different scenario, including 
other synchrotron or inverse Compton components and 
their spectral evolution with time, can make the spec- 
tral shape of the GRB prompt emission more complex 
in the Fermi LAT energy range as for example recently 
observed in the case of GRB 090902B (Abdo et al. 2009). 
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